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Abstract 

The shear (rj) and bulk (C) viscosities are calculated in a quasiparticle relaxation 
time approximation. The hadron phase is described within the relativistic mean 
field based model with scaled hadron masses and couplings. The quark phase is 
treated in terms of the heavy quark bag model fitted to the lattice data. A two- 
phase model allowing for the first order phase transition from the hadron phase 
to the strongly coupled quark gluon plasma is constructed by means of the Gibbs 
conditions. Temperature and baryon density dependence of the calculated viscosity- 
to-entropy ratios (ry/s, C/s) are analyzed and compared with those obtained in other 
models. Special attention is paid to the behavior of viscosity coefficients near the 
critical temperature, from both hadron and quark-gluon side. Effects of resonance 
widths on viscosities and viscosity-to-entropy ratios are estimated. 



1 Introduction 



The study of the transport properties of nonequilibrium systems not far from 
an equilibrium state has a very long story. Methods for the calculation of 
transport coefficients were probed, e.g., in description of nonrelativistic classi- 
cal gases [Ij , liquids and glasses [2] , relativistic gases [3] , cold atomic gases [1] , 
Fermi [5] and Bose [6] liquids. In the past, transport coefficients for the nu- 
clear matter were also studied [TllHllQllTOlITT] . A knowledge of various transport 
coefficients is required also in astrophysical problems such as the entropy pro- 
duction in the universe, the electro- weak baryogenesis [12], for the description 
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of various phenomena in supernovas and neutron stars [TSj. A recent paper [T3] 
considered shear viscosity effects in excited atomic nuclei. With accessibihty 
of heavy-ion colhsions the possibihty for the creation of a new state of matter, 
the quark-gluon plasma (QGP), has been offered. In this respect, there ap- 
peared growing interest in the calculation of the transport coefficients of the 
QGP [^15|16|17|18|19j . Transport coefficients in hot gauge theories were con- 
sidered in Refs. |20|21|22|23|24|25|26|27j . Attention has mostly been focused 
on ultrarelativistic theories, where the scale is set exclusively by the tempera- 
ture T. These results are valid for T much higher than the critical temperature 
Tc of the deconfinement phase transition. Near the critical temperature the 
highly nontrivial effects of the strong coupling should be important. 

Recently, the interest in the transport coefficient issue has sharply been in- 
creased in heavy-ion collision physics. Large values of the elliptic flow V2 were 
observed at the Relativistic Heavy Ion Collider [25]. This finding indicates 
that the created QGP behaves as a nearly perfect fluid with a small value of 
the shear viscosity-to-entropy density ratio, ry/s, which was confirmed by non- 
ideal hydrodynamic analysis of these data [29] . It was claimed [2T|27f30] that a 
new state produced at high temperatures is most likely not a weakly interact- 
ing QGP, as it was originally assumed, but a strongly interacting quark-gluon 
plasma (sQGP). This interest was also supported by a new theoretical per- 
spective, namely. A/" = 4 supersymmetric Yang-Mills gauge theory using the 
Anti de-Sitter space/Conformal Field Theory (AdS/CFT) duality conjecture. 
Calculations in this strongly coupled theory closely related to QCD give for 
the Tj/s ratio [2T|25f31] 



This result was confirmed in subsequent investigations [32|33f34j . After evalu- 
ating ?7/s in several strongly coupled theories it was conjectured that 'f]/ s = ^ 
is in fact a lower bound on the ?7/s-ratio in all systems. For certain materi- 
als, e.g. helium, nitrogen and water, the ?7/s-ratio has a minimum at the 
phase transition [35]. Thus, there appeared belief that QCD nearly saturates 
mentioned minimum near the critical point of the phase transition [36|37] . 
These expectations agree with estimates of the shear viscosity obtained within 
the lattice QCD [55|39f4U|l^ . Actually, there is an extra contribution of soft 
modes in the vicinity of the phase transition critical point which may undergo 
a weak divergence, e.g. for the first order liquid-gas phase transition {H class 
of universality) one expects the behavior 77 oc |T — Tc\~'^^^'^ ~ |T — Tc|~'''°^^, 
see [12]. However, such a weak divergence can manifest itself only in a very 
narrow vicinity of the critical point. 

The bulk viscosity is much less studied than the shear viscosity. At high tem- 
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peratures, when coupling is weak and theory is nearly conformal, the bulk 
viscosity is expected to be very small [T7|43f44f45f46] . In liquids the shear 
and bulk viscosities are usually of the same order of magnitude. However, in 
some cases the bulk viscosity can be significantly higher than the shear vis- 
cosity. This is true in the presence of soft slowly relaxing collective modes. For 
example, such behavior occurs in the vicinity of the phase transition critical 
point or at the crossover |46|47|48|49|50|51] . In this case divergence of the 
bulk viscosity is rather strong, ^ oc |T — Tc\~'^^~^" ~ |T — Tc\~^'^ (the critical 
index z is known from the if-model and u, a from 3d Ising model universal- 
ity class) [12]. Lattice calculations available for the gluon plasma are not in 
disagreement with the expectation of an increase of the bulk viscosity to the 
entropy density ratio, (/s, toward the QCD phase transition critical point 
from above but error bars are still very large |10|52] . There are arguments, see 
[SniSl], that the dynamics of the first-order phase transition is controlled by 
finite values of the kinetic coefficients. 



In the modeling of the strongly interacting matter, interactions are often 
treated within the quasiparticle approximation. Quasiparticle models [55|56] 
describe the lattice data rather appropriately above the critical temperature 
Tc. Relativistic mean-field based quasiparticle models were successfully applied 
to describe the hadron phase [FTII^ . In Refs. the shear and bulk vis- 

cosities of the quark phase were calculated within the quasiparticle approach 
in the relaxation time approximation in the case where the effective masses of 
the constituents depend on the temperature and on the baryon density. 



In this paper, we extend the investigation of shear and bulk viscosities within 
the quasiparticle models in the relaxation time approximation [T7|45f59] . In 
Sect. [21 we describe the hadron phase (T < Tc) in terms of the quasiparticle 
relativistic mean-field-based model with the scaling hadron masses and cou- 
plings (SHMC) [STfSH] . For the QGP phase in Sect. |3]we use the "heavy quark 
bag" (HQB) model which rather appropriately fits the lattice data. The equa- 
tion of state for the two-phase SHMC-HQB model is constructed in Sect. Si 
In Sects. [5] and [6] we evaluate the shear and bulk viscosities for quasiparticle 
collisions in the relaxation time approximation and then present results of nu- 
merical calculations first for the hadron and then for quark-gluon phases and 
compare them with previously obtained ones. In Sect. [TJ we estimate effects of 
finite mass- widths of resonances on the viscosities and viscosity-to-entropy ra- 
tios. The conclusion remarks are given in Sect. |Hi Some details of calculations 
are deferred to the Appendices A-C. 
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2 The SHMC model setup 



Consider hadronic matter in thermal equilibrium. Within our relativistic mean- 
field SHMC model we present the Lagrangian density of the hadronic matter 
as a sum of several terms: 



^ — C,ha.r + -^^MF + -^^cx • (2) 



The Lagrangian density of the baryon component interacting via cr, p mean 
fields is as follows: 



tiG{bar} 



■ (3) 



The considered baryon set is {b} = N{93S), A(1232), A(1116), S(1193), 
S(1318), E*(1385), S*(1530), and (](1672), including antiparticles. The used 
(T-field dependent effective masses of baryons are [57ll58f60j 

ml/rrib = ^biXa(^) = ^ - Qab Xa (r/rrib , b E {b} . (4) 

In Eqs. ([3]), (Hj) g„b, g^b, gpb are coupling constants and Xa{o-), Xtoicr), Xp{<^) 
are coupling scaling functions. 

The a-, u-, p-meson contribution is 
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' - d^Pf, . 

The mass terms of the mean fields are 

'^mMm = |$m(X<TO-)| , {m} = (T, W, p . (6) 

The dimensionless scaling functions $6 and as well as the coupling scaling 
functions Xm depend on the scalar field in the combination Xo-(c") c- Follow- 
ing [60] we assume approximate validity of the Brown-Rho scaling ansatz in 
the simplest form 
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$ = $jv = $^ = $^ = $^ = 1 - /, f = g^N Xa (r/rriN (7) 

with Xa = ^a- 

We keep a standard form for the potential of the non-hnear self-interaction U 
used in relativistic mean-field models, now expressed in terms of /-variable: 



The third term in the Lagrangian density ([2]) includes meson excitations 

■^ex ^ ^ -^bos) (9) 

bosSjex} 

{ex} = vr; K, K; r/(547); a', uj', p; K\ i?*(892), ^(958), 0(1020). 

The knowledge of the Lagrangian density ([2]) defines unambiguously the energy- 
momentum tensor T^^ (Greek indices n,v run 0,1,2,3). The energy density 
E and pressure P are given by the diagonal terms of the energy-momentum 
tensor E = (Too) ' = | {Ta), Latin index i = 1,2,3. The pressure can be 
presented as a sum of the mean-field terms and contributions of the baryons 
and meson excitations [FTfSS] : 

P[f,uo]=PMF[f,u;o]+ E Pb[f,u;o]+ E nos[/,cuo] • (10) 

bS{bar} bosG{cx} 

Below we consider isospin-symmetric nuclear matter (A^ = Z). Then there are 
only a and mean field solutions of equations of motion. Therefore further 
we set the mean fields cJ = and = 0. The value of the mean field uJo{f) 
is found by minimization of the pressure and then it is plugged back in the 
pressure functional. The latter becomes a function only of /. The equilibrium 
value of / can be found by subsequent minimization of the resulting pressure 
in this field. In a self-consistent treatment [58], equations of motion for the 
mean fields render 



^ P[f, Uo] = , ^ P[f, uoif)] = ^ P[f, cooif)] = . (11) 



duo df df 

Since the boson excitation term depends on the mean fields, its minimization 
produces extra terms in the equations of motion for the mean fields. This self- 
consistency of the scheme allows us to be sure of thermodynamic consistency 
of the model. 



5 



It is convenient to introduce renormalized constants 



= "^^^^^ (12) 
mm 

and, instead of Xm, consider other functions 

Vmif) = ^lif)/Xl{f) for{m}=u;,p, r/.(/) = 1. (13) 



In terms of this new notation the contribution of mean fields to the pressure 
am is : 



PuAf, c^o] = VM) - UU) + X. ^of . (14) 



The density of particle species a is given by 



7 / PCX] 

n^iT, fia) = ^ d\p\ F:^{k, T, /ij . (15) 

/TT^ JO 

The corresponding contribution to the pressure in (|TOl) is as follows: 



P.[/,^o] = yy^^^. (16) 

" 

Here a = (6 G {bar}, bos G {ex}), "6" is the baryon and antibaryon, "bos" is 
the boson and antiboson excitation as described above; Va is the degeneracy 
factor; 



-1 



;i7) 



is the quasiparticle occupation for fermions (+) and bosons (— ) with the 
quasiparticle energy 



effective mass m^^^^* and the gauge-shifted values of chemical potentials 
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,,* +ciar,, , j.str ,, j.vec ^ , , ,, j.vec ^* , , 

A'-bos = ^bos A'-str " ^bos S'w.bos '^0 = A'-bos " ^bos fl'w,bos '^O- (19) 

Baryon (antibaryon) quantum numbers t}^*^'", tf'', t™'^ are the baryon charge, 
strangeness and vector charge, t^os is the vector charge of the boson (antibo- 
son) excitation; gf* are effective couphng constants. We use the same values 
of the effective couphng constants and masses as in [58] . One should speciahy 
note that as in Ref. [58] here and below the coupling constants gabi except for 
nucleons, are additionally suppressed by factor 1/10 compared to values which 
are usually used, see Fig. 3 in [58]. This allows us to fit the lattice results for 
the pressure of the quark phase up to T ~ (220 ^ 230) MeV with our SHMC 
model. 

Summing up over all relevant hadron species we obtain the total baryon charge 
and strangeness in the hadron phase 

ha na{T,fIa) , (20) 

a£{h} 

nstr(T,/ibar,/^str) = Y Sana{T,Ha), (21) 
ae{h} 

where the sum is taken over all hadrons {h} = {b} + {bos}. 

The pressure of boson excitation^l] is the pressure of the ideal Bose-gas of 
quasiparticles in the mean fields with effective masses for a-, u-, p- field exci- 
tations and K mesons, see (fT6|) - f lT9|) . For other particles we use bare masses. 

In order to get pp^'"*, we should expand total pressure P[a,ujo{a)] in a' = 
a — 0"'^'. The term linear in a does not contribute due to the subsequent 
requirement of the pressure minimum in a'^^, dP/df = 0. The quadratic term 
produces effective a - particle mass squared. 



( part*N2 _ _fPMFW^(M^ _ C?^Pmf[/, ^o(/)] / df_Y 

^'^'^ ' ~ da^ ~ dp [daj ■ ^ ' 

Keeping only quadratic terms in all thermodynamical quantities in fluctuating 
fields we disregard the boson excitation and the baryon contributions in (!22|) . 
Within our approximation the effective masses of u- and p- excitations prove 
to be the same as those following from the mean-field mass terms 



^ We pay attention to a misprint in Eq. (33) of [58] for the kaon contribution. In 
denominator enter ojk[^o = -Ro = 0]. 
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m 



part* 



part* 



mp\^M)\ ■ 



(23) 



We compare the results of our SHMC model with those for the ideal gas (IG) 
of free particles. In the IG approximation the density of the particle of the 
species a with the spin-isospin degeneracy factor z/^ is 



is the momentum distribution for fermions (upper sign +) and bosons (lower 
sign — ) in the IG limit. Here the chemical potential 



is related to the baryon //bar and strange /igtr chemical potentials which control 
the conservation of the baryon and strangeness charges, respectively; 6^, Sa are 
the baryon/antibaryon charge and strangeness for the given species. 

In Fig. [1] (left panel) we show the square of the speed of sound = dP/ de at 
/^bar = for the SHMC model (solid line) and compare the result with that for 
the IG model with different species number: with the same hadron set as in 
the SHMC model (long-dashed line), for a vr + p mixture (dash-double-dotted) 
and for a purely pion system (dash-dotted). Short-dashed line shows for 
the Hagedorn model [i.e. an IG of hadron resonances whose mass spectrum 
is assumed to have the Hagedorn form p{m) = m~" exp(m/TH) where a = 4 
and Th is interpreted as an upper bound of the hadron temperature [HI]E2])- 

As is seen, for the pure pion IG the monotonously increases with increase of 
the temperature approaching the ultrarelativistic limit = 1/3 at high tem- 
peratures. For the pion-rho meson mixture, the cj. exhibits a shallow minimum 
at T ~ 170 MeV. The minimum (in the same temperature region) is getting 
more pronounced for multi-component systems (see long-dashed curve). At 
T < 50 MeV the pion contribution is a dominant one, thereby all curves co- 
incideH. The curves for the SHMC model and IG calculated with the same 
hadron set coincide for T < 100 MeV. At T > 50 MeV heavier mesons start 
to contribute that slows down the growth of the pressure and then results in 

^ Note that within the SHMC model pions are treated as an ideal particle gas 




(24) 



where 




(25) 




(26) 
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Fig. 1. Left panel: The sound speed squared in hadron matter as a function of the 
temperature at zero baryon chemical potential. Solid line ~ calculation within the 
SHMC model. Long-dashed line - calculation for the ideal gas of particles performed 
with the same hadron set. Dash-dotted and dash-double-dotted lines are results for 
the purely pion ideal gas and ideal pion-rho gas, respectively. Calculation results 
within the Hagedorn gas model [61162] are plotted by the short-dashed line. Right 
panel: Temperature dependence of effective masses for nucleon, w- and p-excitations 
(solid line) and the c-meson excitation (dashed line) within the SHMC model for 
two values of the baryon density. 

decrease of c^, contrary to the case of the one-component pion gas. At still 
higher temperatures heavier species start to contribute in the case of multi- 
component systems resulting in a minimum in c^(T) for IG (see also [62]) and 
SHMC models. In the temperature range 50 MeV < T < 150 MeV the decrease 
rate of cj. is slower in the Hagedorn model than in the SHMC and IG models. 
Only for T ^ Th > 150 MeV a Hagedorn regime is reached, in which more 
and more energy goes into forming massive resonances and therefore, finally 
Cg drops to zero. As mentioned above, the use of the suppressed coupling con- 
stants Qa^b, except for nucleons, guarantees that at < T < (220 -r- 230) MeV 
the SHMC model results for the pressure are in agreement with the lattice 
results [SB]. Within the SHMC model with this modification, gets a deep 
and rather narrow minimum at a critical point T ^ ~ 180 MeV caused 
by a sharp decrease of the in-medium hadron masses at these temperatures 



^ The minimum of the speed of sound (at T = Tc — 180 MeV) can be associated 
with a kind of phase transition, e.g. with the hadron-QGP crossover, as it follows 
from the detailed analysis of the lattice data, see [63] . 
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(see the right panel in Fig. [H where effective masses of the nucleon, u-, p- 
and cr-excitations are presented). Note that in the m-truncated Hagedorn-gas 
model, remains finite at T = ~ T^: For m < M = (1.5 -7- 2.5) GeV the 
results [22] become very close to those for the IG model presented in Fig. [T] 
The specific behavior — )■ arises in the Hagedorn model only for m — 00. 



3 The HQB model 



For the QGP phase of the IG of the massive quarks, antiquarks and gluons 
{q} = q, q, g, see [56] we define the density of conserving charges similarly to 
Eqs. dni): 



n 



HQB 
bar 



(T,/ibar,/istr) = J2 ^an^a^T.Ha) , (27) 



aG{g} 



n. 



HQB 
str 



(T,/ibar,/^str) = ^an^a^T.Ha) ■ (2J 



ae{g} 



With nonperturbative effects associated with the deconfinement transition 
included into a constant vacuum pressure 5, the total energy density and the 
pressure become 



£«Q^(T, /ibar, /istr) = E ^f^^, /i J + 5 , (29) 

P^^^{T, /ibar, /istr) = E P'^^T, fi^) - B , (30) 
aG{9} 

where for gluons (a = g) Ug = 16. The gluon chemical potential Pg = 0; 
thereby, gluonic energy density e and pressure P depend only on temperature. 

Effective masses of quarks and gluons are generally assumed to be temperature 
and density independent and are treated here as free parameters. From the 
fit of thermodynamic quantities to available lattice data for the Nf = 2 + 1 
system [5^64f65j we get quark masses m„ = nid = 100 MeV, = 450 MeV 
and rUg = 600 MeV at B = (215 MeV)^ 
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4 Equation of state in the two-phase model 



Following a common strategy of the two-phase model [66], we determine here 
the deconfinement phase transition by matching the EoS of the SHMC model 
and the HQB model for quarks and gluons. 

The phase equilibrium between the QGP and the hadronic phase at the 
first order phase transition is determined by the Gibbs conditions for ther- 
mal (T^QB = 2-SHMC)^ mechanical {P^^^ = pSHMC) chemical = 
A^bar^*"' y^st?^ = A^st™*") equilibrium. The chemical equilibrium between differ- 
ent components is automatically satisfied owing to representation fl2^ for the 
chemical potential. At given temperature T and baryon chemical potential 
/ibar the strangeness chemical potential /Xstr is determined by fixing the net 
strangeness of the system to zero. So, in the Gibbs mixed phase at the fixed 
total baryon charge density nbar the following set of equations is to be solved : 



pSHMC^j^^ /ibar; /^str) — P^^^(T, /ibar; /^str) ; 

nbar(T,/ibar;/istr) = « nba^^ (T, /ibar; /^str) + (!-«) Tib^^l^; /^bar; /^str) ; 
= a nf^^{T, /ibar; /istr) + (!-«) n^^^^{T, /ibar; /istr) ; (31) 

where a = V^^^/V is a volume fraction occupied by plasma, i.e. the coexis- 
tence phase is not a homogeneous one. The boundaries of this mixed phase are 
found by putting a = 1 {the hadron phase boundary) and a = [the plasma 
phase boundary). Eqs. (13T|) should be solved at every phase point (T, rzbar) 
of the coexistence region. It results in that the quark-hadron boundary of a 
system conserving baryon charge and strangeness becomes a critical surface 
(not a line in the T — /ibar plane as for a single charge conservation, but some 
stretched-out area p7].) 

Predictions of the hadronic SHMC model with suppressed couplings agree 
with the recent lattice data for the reduced energy and the pressure up to 
about T ~ 220 ^ 230 MeV (see Fig. El left and right). As follows from Fig. El 
at higher temperatures quark-gluon degrees of freedom should be taken into 
account. Thereby, further we use the SHMC-HQB mixed phase model. In the 
/ibar = case the energy density of the SHMC-HQB mixed phase model suffers 
a jump at the temperature T = Tc ~ 180 MeV (the left panel in Fig. E|)- 
This value of the critical temperature is consistent with that calculated on 
the lattice at rzbar = [6lf68] . However the lattice data demonstrate a sharp 
crossover transition at this point. We found a jump in the latent heat about 
Ae ~ 0.85 GeV/fm'^ (Ae/T^ ~ 6). The presence of the jump in e at 
indicates the first order phase transition. From the right panel of Fig. El we see 
that the pressure undergoes a continuous monotonous increase in temperature 
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Fig. 2. Reduced energy density and pressure at ^bar = obtained in the two-phase 
SHMC-HQB model. The energy density jump at Tc ~ 180 MeV marked by the 
dotted line corresponds to the mixed phase. The hadronic SHMC results above 
Tc are plotted by dashed lines. QCD lattice results with the time extent Nt =8 
for the p4 and asqtad actions are plotted by empty diamonds and filled circles, 
respectively [68] . 

in agreement with the construction of the mixed phase. 

In Fig. [3] we demonstrate the reduced entropy density as a function of temper- 
ature. The lattice data using two types of the action are plotted by points for 
Nf =8. The results of the two-phase model and purely hadronic SHMC model 
are shown by solid and dashed (prolonged for T > T^) lines, respectively. The 
SHMC model with suppressed couplings continues to be in close agreement 
with the lattice data till T ~230 MeV. The HQB model agrees reasonably 
with the lattice data [68] for all temperatures above Tc except a narrow vicin- 
ity of the critical point where in the lattice QCD we deal with the crossover 

at /ibar =0. 

Reasonable agreement between the SHMC-HQB model and lattice data is 
observed not only for the reduced energy, pressure and entropy but also for 
the trace anomaly (e — 3P)/T'^, for both zero and finite values of the baryon 
chemical potentials, as demonstrated in the left and right panels of Fig. |H As 
is seen, only at large temperatures {e — 3P)/T^ — )■ 0, i.e. the system becomes 
a conformal one. The conformal invariance in QCD is essentially broken at 
T < 2Tc. Therefore, we may expect here differences in predictions of conformal 
theories for different quantities, e.g. for rj/s, at T < 2Tc. We stress that the 
conformal regime is not reproduced by the purely hadronic SHMC model, 
valid roughly for T < 230 MeV. 

Values of the ratio P/e and the speed of sound Cg essentially determine the 
evolution of a system, and as it will be shown below, the last quantity is 
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Fig. 3. Reduced entropy density at /^bar = obtained in the two-phase SHMC-HQB 
model. The entropy density jump at Tc ~ 180 MeV marked by the dotted line 
corresponds to the mixed phase. Notation is the same as in Fig. [2j 




Fig. 4. Trace anomaly at //bar = and /Xbar = 530 MeV. The jump is marked by 
the dotted line. The hadronic SHMC results are plotted above Tc by dashed lines. 
Lattice QCD data are from [69] (filled squares), and from [68] (circles and diamonds, 
notation is as in Fig. [2]). 

an important ingredient entering into the expression for the bulk viscosity 
coefficient. 
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Fig. 5. The P/e ratio demonstrating the presence of the softest point and the speed 
of sound squared in the baryonless matter are presented. SoUd hnes in both panels 
are results of our two-phase model calculations. Dash-dotted lines continuing the 
solid lines at T > Tc demonstrate calculation results of the SHMC model. Left: 
empty diamonds and filled circles are the Nf =2+1 lattice data for the P/e ratio 
calculated with the p4 and asqtad actions, respectively [68]. The dashed line is 
the high-T approximation (I32p of lattice data with improved staggered fermions 
and almost physical masses for Nf = 2 + 1 [6lj- Right: Diamonds and squares are 
lattice QCD data for Nf =2 with Wilson fermions [63] and Nf =2+1 with staggered 
fermions [70], respectively. 



In Fig. [5] the energy density dependence of the P/ e ratio for the given two- 
phase SHMC-HQB model EoS is presented alongside with the temperature 
dependence of c^. Both quantities {P/e and c^) have close physical meaning 
and coincide in the ultrarelativistic limit = P/e = 1/3. This limiting trend 
is observed in Fig. [5] (left and right panels) at high temperatures and energies. 
In the high energy density region the SHMC-HQB results agree reasonably 
well with the lattice data. Though the cited above lattice results were ob- 
tained in different discretization schemes for 2- and (2+1) flavor QCD using 
different quark masses and lattice spacing, thermodynamics in the high tem- 
perature phase is not strongly influenced by discretization errors and is rather 
insensitive to changes of the quark mass. Nevertheless, one should be cautious 
treating these lattice QCD results for not too high temperatures, especially 
for T < Tc- Thereby, our point here is to use the HQB model for T > T^. 
that fits the lattice results well for T > 2Tc and to apply our SHMC model 
in the hadronic sector at T < Tc, instead of using the corresponding lattice 
values for T < Tc. It is of interest to note that the SHMC model reasonably 
reproduces both P/e and at the temperature above till about 230 MeV, 
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in agreement with results presented in Figs. [2|3|^ 



Lattice data [M] on the P/ e ratio obtained for the pion mass ttlt^ ^ 220 MeV 
are well approximated in the high temperature/energy range 1.3 < e^^^ < 6 
GeV/fm^ as 

^ = ifo.964 ^i^T^V (32) 

£ 3 V 1 + 0.26 e W/GeV / ^ ' 

This approximating curve is plotted by the dashed line in Fig. [5] (left panel). 
We see that different lattice data presented in this figure agree rather well 
with each other. Since = dP/de, the approximation fl32l) can also be used 
to reproduce the speed of sound. It is seen that in contrast with a single 
component hadronic gas (cf. Fig. [T] and Fig. both and P/e exhibit a 
minimum near the phase transition. In other words, the EoS in our two-phase 
SHMC-HQB model is getting softer near Tc and the system evolves slower at 
the "softest point", i.e. at the minimum of the P/e ratio. As was noted above, 
for hadronic sector T < the lattice results predict much lower values of the 
P/e ratio, as compared to our SHMC model. Being in discrepancy with lattice 
results for the low-temperature hadronic component, results of our two-phase 
and SHMC models agree reasonably well with other hadronic models, as it is 
demonstrated in Fig. [TJ 

The following remarks are in order. As was found within the standard RMF 
model [7T] the specific heat exhibits a sign of critical behavior. The SHMC 
model with the standard choice of couplings demonstrates a crossover at a 
pseudocritical temperature T ~ Tc for /ibar = in agreement with the lattice 
data. In this case we observe a sharp peak in the heat capacity which is typi- 
cal for the strong crossover transition [57] . In our case when couplings are 
suppressed, the peak is smoothed that reminds about a weak crossover. Our 
two-phase SHMC-HQB model shows the first order phase transition. This can 
be treated as an advantage of the SHMC model and a disadvantage of the 
two-phase model. However the results of the SHCM model begin to disagree 
with the lattice data for T > 220 — 230 MeV, whereas the two-phase model 
fits well the lattice data in the high temperature region. Thus we can't use the 
SHMC model for high temperatures without its significant modification. We 
refuse of doing such a modification bearing in mind that the high tempera- 
ture regime should be treated in terms of the quark-gluon degrees of freedom. 
The mentioned lack of our two-phase model can be removed by a more careful 
treatment of the near-critical region. On the one hand, in the framework of the 
model with the first-order phase transition our consideration of the vicinity of 
the critical point is obviously oversimplified, since the Gibbs condition's treat- 
ment disregards dynamical effects in the phase transition. The latter effects 
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are important for the description of the vicinity of the critical end point. These 
dynamical effects are manifested at temperature below the critical end point 
in all models with the EoS of the van der Waals type, see [MfTS] . In particular, 
an overcooled state to be possibly formed at the first order phase transition 
may result in a mechanism of fast hadronization [73] . Thus one way to better 
treat the near-critical region is to incorporate non-trivial dynamical effects, 
e.g. in the framework of the van der Waals model EoS. On the other hand, 
one could smooth the thermodynamic quantities in the near-critical region in 
such a way that the EoS would correspond to the crossover for sufficiently 
small values of /ibar, see [71]. We will return to the consideration of the critical 
effects in the subsequent work. 



5 Shear and bulk viscosities in the hadron phase 

5. 1 Collisional viscosity in the SHMC model. Derivation of equations 



Following Sasaki and Redlich [59], we derive expressions for the shear and 
bulk viscosities in the case when the quasiparticle spectrum is given by the 



expression E{p) = yp'^ + m*2(T, /i). We perform a similar derivation, but in 
the presence of mean fields a and uq. In the latter case one should additionally 
take into account that quasiparticle distributions depend on the mean fields. 

We start with the expression for the energy-momentum tensor: 



T''' = T^^+ E E ^b^s, (33) 

fe€:{bar} bos6{cx} 

where the mean-field contribution is as follows 



= {I [mf a' - m*J u'^ + U{a)) (34) 

with m* and given by Eq. (|6]), and where we dropped the terms quadratic 
in gradients. 

The quasiparticle (fermion and boson excitation) contribution is 



= / dV |M_^«Mi7;| . (35) 
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Here 



= Qu^b Xu. w^, X^^ = gl^^^ w^. 

Note that the contribution to the energy-momentum tensor T^'^ is not sym- 
metric with respect to the interchange of indices. The tensor can be sym- 
metrized following a general rule, e.g. see [75]: One may add to T^'^ an extra 
term T'^^, such that d^.T'^'" = and that T'"' + T''"' = T"'' + T'"''. Thus the 
asymmetry of T^'^ does not influence on observables. In addition, to calcu- 
late viscosities we need to deal only with spatial components of the tensor, 
T**^. Since = {Xq, 0, 0, 0) in the mean field approximation used by us, up 
to second gradients we have T^'^ ~ T^*. Therefore we will not perform this 
symmetrization procedure in the given paper. 

The quasiparticle distribution function Ff, for baryon components in the pres- 
ence of mean fields fulfills the Boltzmann kinetic equation [76] , 



{p'td, - g.bP.uj'"^^ + ^d'^ml^^ Fb = StFb- (37) 
with = 5{pI - m*2)F^-^q-(pl„x,). 

The local equilibrium boson or baryon distribution is given as follows: 



where we suppressed u'^ terms for \u\ <^ 10 Here the upper sign (+) is for 
fermions and (— ) is for bosons, and the vector particle charge is t^*^^ = ±1 or 
0; Ouja 7^ only for a G bar and kaons in our model. Considering only slightly 
inhomogeneous solutions and using |n| ^ 1 we may drop the terms oc and 
(X uS/ojq in the kinetic Eq. (137|) . Then kinetic equations for boson and baryon 
components acquire ordinary quasiparticle form 



dFa d Fad Fa d Fad Fa p'i dFa ^ ^ 

+ -T^l^ - 1^1^ = ^1^ = StFa, (39) 



dt dpa dfa dfa Opa Fa Ox^ 



^ For u ^ there appear mean field solutions with uj ^ 0. These however yield 
<C 1 terms. 
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where = {Ea{pa,ra, a,u),pa). We used that dEa/dpa = Pa/Ea- Since cal- 
culating the viscosity, we need only terms with velocity gradients, we further 

put dEjdfa = (9E„/9/iJVa/ia + {dEjdT)WaT = 0. 

In the relaxation time approximation 



StFa = -6FjT^, 5F, = F,-F':--'^-. (40) 

Here denotes the relaxation time of the given species. Generally, it depends 
on the quasiparticle momentum pa- 

The averaged partial relaxation time fa is related to the cross section as 



f,"i(T,/i) = ^n,,(T,^) (tW^L'KaO), (41) 

a 

where n^' is the density of a'-species, = / d cos 9 da {aa — )■ aa')/dcos6 (1 — 
cos 6') is the transport cross section, in general, accounting for in-medium 
effects and v^^' is the relative velocity of two colliding particles a and a in case 
of binary collisions. Angular brackets denote a quantum mechanical statistical 
average over an equilibrated system. In reality, the cross sections entering 
the collision integral and the corresponding relaxation time in (l40l) may 
essentially depend on the particle momentum. Thus, averaged values given 
by Eq. fl41l) yield only a rough estimate for the values t^^ which we actually 
need for calculation of viscosity coefficients, see below Eqs. (H8|) and fHOj) . 

In the relaxation time approximation from Eqs. fl39l) . pUj) we obtain 



nploc.eq. 

and then the nonequilibrium correction to the energy-momentum tensor (1331) 
becomes: 



6T 



E/^rjr.M^^p^S.Fj (43) 



loc.eq. 



5a — — + 5(jJq 



da loc.eq. ^UJ^ 



loc.cq. 



Considering small deviations from the local equilibrium, we may keep in (143 
only first-order derivative quasiparticle terms oc 9^. 
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The shear and bulk viscosities expressed through traceless nondiagonal and di- 
agonal parts of the nonequilibrium correction to the energy-momentum tensor 
cirG clS follows: 



"C %V ■ U - 7] W, 



(44) 



Wki = dkUi + diUk - -5ki diu\ 
As before, here Latin indices run 1, 2, 3. 



To find the shear viscosity, we put i ^ j in f l44j) and use that in this case the 
variation of the second and third terms in ( l43l) yields zero after integration 
over angles. To find the bulk viscosity, we substitute i = j m. (jS]) and use that 
= 3Peq- As foUows from equations of motion ( ITTi) . variation of the second 
and third terms in (143|) again yields zero. We put m = in final expressions 
but retain gradients of the velocity. 



Taking derivatives 9F^°'=-"'^-/(9x^ in Eq. (02]) with the help of Eqs. (1791) - (IHTD 
from Appendix A, we find the variation of the total energy-momentum tensor 
as the function of derivatives of the velocity 



(45) 



with the upper sign (— ) in the blocking factor for fermions and the lower one 
(+) for bosons in accordance with Eqs. f l25|) . fl38l) . 



(46) 



Qa 



X 



dP 



-0 



0^ 



/ dP 

V^'^bar 

djEg + X°) 
dT 



d{Ea + XO) 



±bar 



bar 



str 



— /^b£ 



'dP\ 

de 



d{E, + 



0^ 



Z^str" 



d{Ea + X: 



bar 



str 



(47) 



^ Since we are not interested in calculations of the heat conductivity, we again 
suppress the VT, Vfi terms and also put u = keeping only the corresponding 
derivative terms. 
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Finally, we obtain the shear viscosity 

'^ = I^i:/<'rr„|KMiTFr)l. 

Correspondingly, the bulk viscosity is 



(48) 



(49) 



To get a final expression for the latter quantity, one should use the energy 
conservation obtained from Eq. f l45l) by demanding ^T^*^ = 0, 



E / ^ r [r, {E, + xl) Fr (1 T qa 

+ fa {Ea - X°) Fr (1 T 







(50) 



where the second term corresponds to antiparticles with the relaxation time 
fa- After substitution of Eq. (l50l) in f H9|) . the bulk viscosity can be presented 

as 

- ^ E / dr (-^ + Fr (i T f:'') g., 

where Qa is defined by Eq. ( H71) . 

For X° = 0, A^str = our results coincide with those obtained in [59]: 



(51) 



Ca = y dF (1 T X 



(52) 



X 



3F, 



9P 



9n 



9e J 



bar . 



dE 



ibar 



9/i 



bar 



(53) 



9F„ dEa 

i^n — 1 -::rz=. /iba 



dT 



where the quasiparticle mass m* can be temperature/density dependent. Mak- 
ing substitution of the on-shell mass instead of m* and F^^ = F^^ in Eqs. fl52|) . 
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( !53|) we arrive at expressions for the gas of non-interacting particles. For the 
gas of quasiparticles with values m* depending only on the temperature with 
the help of Eq. f l50p . Eq. f l53p can be rewritten in a simpler form 



C = i/drrF-(lTF-){(i-c?)E-|^}\ (54) 

In a particular case of the one component gas of material particles in the high 
temperature limit T ^ m for r = f = const viscosity coefficients take the 
simplest form [16]: 



V-^e{T), C^O, (55) 
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where e(T) = ^vr^T^ is the known black-body radiation result for a massless 
Bose fluid. 



5.2 Collisional viscosity in hadronic haryonless matter 



In the relaxation time approximation both shear and bulk viscosities for a 
component "a" depend on its relaxation (collisional) time r^, which should 
be parameterized or calculated independently. Therefore to diminish this un- 
certainty it is legitimate at flrst to flnd the reduced kinetic coefficients (per 
unit relaxation time, assuming r = const, i.e. r = f). 

In Fig. [6] we demonstrate results of various calculations for the reduced shear 
(left panel) and bulk (right panel) viscosities at /Xbar = additionally scaled 
by the \/T^ factor. As we see from the figure, the reduced scaled shear viscos- 
ity of the massive pion gas (dashed line) becomes approximately constant for 
T > 100 MeV. Naturally, this result is close to that obtained in the Gavin ap- 
proximation [16] where a numerical interpolation between the m/T <^ 1 and 
m/T 3> 1 cases was used (dash-double-dotted line in Fig. [6]). The scaling 
is violated for the n — p gas in the temperature interval under consideration 
because the p mass is not negligible even at T ~ 200 MeV. For ( the approx- 
imate l/T'^ scaling property holds for the massive pion-rho gas at T > 150 
MeV. Note that ( =0 for the gas of free massless pions since = 1/3 in this 
case. For the massive pion gas C/T'^ gets maximum for T ~ 60 MeV and then 
begins to decrease for T > 60 MeV reaching zero at large T similar to the case 
of massless gas. The reduced shear and bulk viscosities of a multicomponent 
system calculated in our SHMC model (solid lines) and in the IG model with 
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Fig. 6. The reduced scaled shear (left panel) and bulk (right panel) viscosities 
per unit relaxation time calculated within the SHMC model for a system with 
Mbar =0 (solid lines). The results are compared with those for the massive pion 
gas (dashed lines), tt — p mixture (short dashed line) calculations and with some 
interpolation from massive to massless pion gas (the Gavin approximation |16j . 
dash-double-dotted line) as well as for the IG model with the same large set of 
species as in the SHMC model (open dots). 

the same hadron set (open dots) do not fulfill the scaling law. These mod- 
els include a large set of hadrons. Due to that with the temperature increase 
the reduced shear and bulk viscosities become significantly higher than those 
for the pion gas and the pion-rho gas models. An additional increase of the 
reduced viscosity within the SHMC model originates from significant mass 
decrease at temperatures near the critical temperature. The bulk viscosity of 
a single-component pion system drops to zero both at low and high temper- 
atures and in the whole temperature interval ( « rj, that is frequently used 
as an argument for neglecting the bulk viscosity effects. However, the state- 
ment does not hold anymore for mixture of many species. For example, at 
T ~150 MeV the t]/( ratio is only about 3 in the case of the IG and SHMC 
models. Thus the bulk viscosity effects can play a role in the description of 
the hadronic stage at high collision energies, like at RHIC. Moreover, the bulk 
viscosity can be responsible for such important effect as flow anisotropy. 

For further evaluation of the absolute values of the transport coefficients r]a 
and (a we need to estimate the relaxation time of hadronic species. For that 
we use Eq. (HT]) . We apply free cross sections in the case of the IG based 
model, similar to procedure performed in Ref. [77]. For the SHMC model, 
the in-medium modification of cross sections is incorporated by a shift of a 
"pole" of the collision energy by the mass difference — m* according to 
prescription of Ref. [78]. Due to a lack of microscopic calculations this is the 
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only modification, which we do here. 



The relaxation time for the tttt collisions in the purely pion gas is plotted in 
Fig. [71 The same cross sections as in Ref. [77] are taken since pions within our 
SHMC model are assumed to have free dispersion law, as well as in Ref. [77] . 
Therefore our results (solid and dash curves) coincide with those in Ref. [77] 
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Fig. 7. The pion relaxation time for the tttt system as a function of temperature. 
Solid and dashed lines are our calculation results for quantum Bose and classi- 
cal Boltzmann statistics, respectively. The corresponding points (empty and black 
squares) are taken from the review-article of Prakash et al. [77]. 

(empty and black squares). As is seen, there is only a small influence of quan- 
tum statistics on the pion relaxation time. 

The temperature dependence of the ratio of transport coefficients to the en- 
tropy density is analyzed in Fig. |S] for the /^bar =0 system. Here we make 
emphasis on the high-T behavior to clarify the question how hadronic models 
describe approaching the lower bound of the rj/s = l/Air ratio predicted by 
strongly coupled theories [2T|25|31j . As is seen from Fig. |8] (left panel), the 
SHMC model (solid lines) as well as the excluded volume one (short dashed 
and dash-dotted lines) [75] predict a monotonous decrease of the specific shear 
viscosity r]/s. However, only in the SHMC model the lower bound rj/s ~ I/Att 
is reached at the critical temperature Tc ~ 180 MeV. The hadron model [80] 
describes all the known particles and resonances with masses rrii < 2 GeV in 
terms of an excluded volume model with r = 0.5 fm and additionally includes 
an exponentially increasing number of Hagedorn states for rrii > 2 GeV. It 
is natural that for T <140 MeV the resonance gas model with the Hagedorn 
states reproduces the behavior oi rj/s for the hadron excluded volume model 



23 




Fig. 8. Specific shear (left panel) and bulk (right panel) viscosities {rj/s and C/s) 
of hadrons as a function of T for //bar =0- Our SHMC model results are shown by 
solid lines. The results of the excluded volume hadron gas model [79] are shown for 
two values of hard-core radii r (short dashed and dash-dotted lines). The results 
for the resonance gas model including Hagedorn states |80j are shown by the long- 
dashed lines. The dotted horizontal line corresponds to the lower AdS/CFT bound, 
r]/s = l/47r [H]. 

for r = 0.5 fm. At T ~ results for the shear viscosity calculated in the line 
of classical non-relativistic approximation [80] are consistent with the lower 
bound of the rj/ s ratio. Note that our SHMC relativistic mean-field model does 
not need extra assumptions on the particular Hagedorn states to describe the 
dip in Cg and rj/s ~ l/47r near the critical temperature. 

The specific bulk viscosity C,/ s (see Fig. [HI right panel) behaves quite differently 
in our SHMC and the Hagedorn hadronic models since completely different 
assumptions were used for calculation of the bulk viscosity. Namely, our results 
are based on the quasiparticle coUisional relaxation time approximation while 
in [80], following the QCD sum rules in [50|5T] . the bulk viscosity C, was 
related with the trace anomaly < 6 >t=< >t= £ — 3P: 



C(T) 



9ujo{T) 



T^A<^>^-<^>0 + 16|,,| 



dT 



rp4 



(56) 



where \eo\ =< 9 >o /4 in vacuum. To derive Eq. (1561) . a particular ansatz 
was used for the spectral function at zero momentum. However, recently this 
ansatz has been criticized in [1^52ll81ll5?| . Moreover, Eq. obviously is not 
applicable for temperatures far from T^. 
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Viscosities in the interacting hadronic phase were calculated in Ref. [83] for 
T < Tc when the dominant configuration of QCD with two fiavors of massless 
quarks is the interacting gas of chiral pions. In the large limit for massless 
chiral pions it was obtained that 



P. 

e' 



— c„ 



(57) 



As was first shown by Gavin [16] in pure massless ideal pion gas C = 0. However 
if pions are coupled with themselves or with the matter of other species e ^ 3P 
and Ctt 7^ 0. Eq. (IFTI) is similar to a simplified expression (/rj ~ 15 (| — c^Y 
which is obtained for a photon gas coupled to hot matter [HI] and is also 
parameterically correct for perturbative QCD [35]. This is because 2 — )■ 2 
scattering is the dominant process in both t] and ( calculations. It is however 
not the case, e.g. in the X(j)^ model, where t] is dominated by 2 <H- 4 processes 
while C by 2 — )■ 2. For a strongly coupled N = 2* gauge theory using AdS/CFT 



the scaling is (/rj oc (| — c^)^ 
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Fig. 9. The ratio of the bulk-to-shear viscosity for baryonless hadronic matter calcu- 
lated within the SHMC model (solid line) and for the massive ideal pion gas (dashed 
line). For comparison T-dependence of the ratio calculated following Eq. ()57p (i.e. 
as for chiral pions, but with thermodynamical quantities as for the massive ideal 
pion gas and normalized to the value for the massive pion gas at large T), is shown 
by the dash-dotted line. 

In Fig. [9] we show the ratio (/rj calculated for baryon-less matter in SHMC 
model (solid line) and those computed for the massive pion gas with a free 
dispersion law (dashed). For comparison the (/rj-Tatio calculated according 
to Eq. ( 1571) (i.e., as for chiral pions) but with thermodynamic quantities the 
same as for the massive pion gas is shown by the dash-dotted line. At T > 
pions become almost chiral ones and the T dependence of the pion gas 



25 



indeed can be described by Eq. (1570 . This coincidence of two curves at large 
T allowed us to estimate the value of the numerical pre-factor in Eq. (1571) . For 
Z^bar = the C,/ri ratio in our SHMC model coincides at T < 50 MeV with 
that for the free massive pion gas. It demonstrates that the main contribution 
to viscosities is given here by the pion mode. In the purely pion case the ratio 
C/?7 monotonously decreases with increase of the temperature, whereas within 
our SHMC model it gets a minimum at T ~ 50 MeV and then a maximum 
for T ~ 120 MeV. 

5.3 Collisional viscosity in baryon enriched hadronic matter 

Similarly to the /ibar =0 case, to avoid influence of uncertainties in the value 
of the relaxation time we study at first reduced kinetic coefficients. Transport 
coefficients for a nucleon-antinucleon mixture at finite baryon density were cal- 
culated earlier by Hakim and Mornas within the standard Walecka model [11] . 
Comparison between their (dotted lines) and our (solid lines) SHMC model 
results is given in FigHOlfor two values of the density. Since Ref. [TT] presented 
results for one component (neutron) matter, we performed calculations in the 
SHMC model also for neutron-antineutron matter. Because the p mean field 
is absent in the standard Walecka model, we suppressed the p-term also in the 
SHMC model calculation. We see that the reduced T^-scaled shear viscosity 
behaves similarly in both models but absolute values differ due to different 
effective masses. Only at high temperatures and for ribar — '^o ('^o is the nu- 
clear saturation density) the values of reduced shear viscosities become close 
in both models. Differences in the values of the reduced bulk viscosity calcu- 
lated in two models are larger than for the reduced shear viscosity because C, 
depends on the specific behavior of thermodynamical quantities (cf. Eqs. (1521) 
and (135]) ). In particular, it concerns the high-density region (see the curves 
for TT-bar = 4?7,o), whcu the mass in the standard Walecka model is getting very 
small. 

The temperature dependencies of the reduced T^-scaled shear and bulk vis- 
cosities calculated for the case of the multi-component hadron mixture within 
IG and SHMC models are shown in the left and right panels of Fig. [TT] at 
baryon densities nbar = and Auq. The reduced shear viscosity calculated in 
the SHMC model (solid lines) is close to that in the IG model with the same 
hadron set (dashed lines). Differences in the ?7/(rT^) ratio for the IG and 
SHMC models appear only at high temperatures T > 150 MeV. At T < 100 
MeV the reduced T^-scaled bulk viscosity (right panel) in the IG based model 
proved to be larger than that in the SHMC model. Contrary, for larger T the 
values of the reduced bulk viscosity in the IG model become smaller than those 
in the SHMC model. Differences come from the strong dependence of the bulk 
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Fig. 10. The temperature dependence of the reduced T^-scaled shear (left panel) and 
bulk (right panel) viscosities per unit relaxation time calculated for a one-component 
pure nucleon-antinucleon system at nbar = and 4no. Solid lines are the SHMC 
model calculations. The standard Walecka model results [llj are plotted by dotted 
lines. 
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Fig. 11. The SHMC model predictions of the T^-scaled temperature dependence 
of the reduced shear (left panel) and bulk (right panel) viscosities calculated for 
hadron mixture at ribar = ''^o and Ahq (solid lines). Calculations performed in the 
IG based model with the same hadron set as in the SHMC model are demonstrated 
by dashed lines. 

viscosity ( on the values of thermodynamical quantities (see Eqs. ( 1471) . (149|) ). 
Note that at T >100 MeV and nbar ~ "'^-o the shear and bulk viscosities are 
getting comparable in magnitude. Growth of the relative importance of ( with 
increase of temperature seems to be quite natural because the bulk viscosity 
takes into account momentum dissipation due to inelastic channels and those 
number grows with the temperature increase. 

Now let us proceed to an estimate of the relaxation time for a baryon enriched 
nuclear system. Important peculiarity of the nucleon contribution to the relax- 
ation time at low temperature is associated with the particular role played by 
the Pauli blocking. It means that appropriate multi-dimensional integration 
should be carried out quite accurately with using quantum statistical distri- 
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Fig. 12. The temperature dependence of the partial nucleon-nucleon relaxation 
time for nbar = and Ino- SHMC model calculations with the relativistic Boltz- 
mann distribution function performed following (j4ip with the experimental free 
nucleon-nucleon cross section are shown by solid lines and those with the effective 
cross-section a^N = Ce// =40 mb, by dotted lines, respectively. The full calculation 
for nucleon component at n = no performed with the Fermi distribution function 
is presented by stars. The ideal Boltzmann nucleon-antinucleon gas model results 
are shown by dashed curves. The Danielewicz interpolation given by Eq. (j58p of [9] 
is plotted by dash-dotted lines for Ino and 4no. Squares are results of calculations 
[77j for ntiar = performed within the Boltzmann approximation. 

bution functions. Calculations using the kinetic Uehling-Uhlenbeck equations 
for the purely nucleon system (A^ = Z) in the non-relativistic approximation 
were performed in Ref. [9]. For T < 100 MeV an interpolated expression has 
been obtained: 



850 / nba 



r2 V no 



1/3 r 
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rV2(i + i60/T2)nba 
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Here t^n is in fm/c, T is in MeV, no = 0.145 fm ^. Thus the relaxation time 
demonstrates well known T"^ behavior for a Fermi liquid at T — )■ 0, cf. [85] . 



In Fig|T2]we compare evaluation fl58|) with the relaxation time, as it follows 
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from the SHMC model for the purely nucleon-antinucleon system with = 
Z . Here we use the Boltzmann approximation and free cross sections just 
corrected by the nucleon mass shift (see solid lines). As is seen from FigJT^ 
due to an account for the Pauli principle the Danielewicz's f^^ given by Eq. 
fl58l) (see dash-dotted line) essentially deviates from the SHMC model result 
(computed here within Boltzmann approximation) for T < 100 MeV. The 
SHMC model results are very close to those for the ideal nucleon-antinucleon 
gas (dashed lines) for T < 150 MeV. A partial account of the Pauli principle 
by assuming constant effective cross section a^jf =40 mb (see dotted line) 
improves SHMC model agreement with the estimates fISS]) at low T. Our full 
calculation of f^j^ for a purely nucleon gas that includes Pauli principle plotted 
by stars demonstrates a reasonable agreement with (!58|) . Therefore to simplify 
calculations below we use Eq. (158!) for the partial nucleon-nucleon relaxation 
time TTVAr, to be valid at low temperatures, smoothly matching it (at T ~ 100 
MeV) with the partial nucleon contribution calculated following Eq. (1411) for 
higher temperatures. In the case nbar =0 the ideal Boltzmann gas model results 
(dash curve) are very close to those of |77j (squares) and agree with the SHMC 
calculations (also performed in the Boltzmann approximation, see solid and 
dotted lines) for T < ISOMeV. For higher T the difference with the IG model 
is due to a sharp change of the effective nucleon mass in the SHMC model at 
large temperatures. 

In subsequent calculations we take into account the whole hadron set involved 
into the SHMC model. The relaxation time for every component, except for 
nucleons, is evaluated according to Eq. (HTi) . 

In Fig{T3] (left) the partial contribution of the nucleon shear viscosity for the 
multi-component system is pictured as a function of temperature at ?7,bar = I'^o 
and 4no. Solid and dashed lines are our results for the SHMC and IG models, 
respectively, with the relaxation time calculated following Eq. (!58l) . The dot- 
dashed lines show the Danielewicz analytical fit [9] 



where r] is given in MeV/(fm2c), T in MeV. For both ribar/^o =1 and 4 the IG 
model results are rather close to this fit, provided fj^^ is calculated accord- 
ing to (l58l) . The shear viscosity calculated under the same conditions within 
the SHMC model is higher than that calculated following Eq. (!59|) and in IG 
model, since the former takes into account decrease of the effective masses 
with the increase of the baryon density. Two short-dashed curves at the bot- 
tom of FigllS] (left) represent 77 computed in the SHMC model within the 
Boltzmann approximation for densities ribar = '^o and 4no. As is expected. 



(irOO/T^) (nbar/no)' + [22/(1 + ■ 10"^)] (nbar/rio) 

+ 5.8ri/V(i + leo/T^) , 



0.70 



(59) 
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Fig. 13. The T-dependence of the partial nucleon contribution to the shear and bulk 
viscosity in the multi-component system for different values of the baryon density 
nbar- Solid and dashed lines are our results, respectively, for the SHMC and IG 
models with the same particle set, with the relaxation time calculated following 
Eq. (j58p . Dot-dashed curves are the Danielewicz interpolation (j59p of results [9]. 
For comparison in left panel two short-dashed lines are shown for the SHMC model 
calculations done within the Boltzmann gas approximation at nbar = 4no (upper 
line in low T region) and no- In right panel bold solid and dashed lines correspond 
to TT-bar = 4no and thin ones, to nbar = n-0- 



in the Boltzmann case at low temperatures the shear viscosity only weakly 
depends on the baryon density (cf. short-dashed curves in Fig JT3l) . This has 
been noted earlier in [8]. Our results demonstrate important role played by 
the Pauli principle at low T and show a strong density dependence. 

The model dependence of results is more pronounced for (. In the right panel of 
Fig{T3]we present the partial contribution of the nucleon bulk viscosity for the 
multicomponent system. We see that in the SHMC model this contribution is 
very small and can be neglected. In the IG model the partial contribution, (n, 
proves to be negative, whereas the total value ( remains, certainly, positive. 

The temperature dependence of the specific rj/s and (/s transport coefficients 
is presented in FigHHfor two values of the baryon density Iuq and Auq. The 
results of the SHMC (solid and dashed lines) and IG (dash-dotted and short- 
dashed lines) models differ not very much for the given value of the baryon 
density, except for the shear viscosity at low temperature due to the Pauli 
exclusion principle. The rj/s for the baryon rich matter achieves the lower 
AdS/CFT bound rj/ s = ^ at smaller T < Tc for the higher ribar- 
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Fig. 14. The T-dependence of the specific viscosity coefficients r]/ s and C,/ s for two 
values of the baryon density nbar- Solid and dashed lines are our results for the 
SHMC model for no and 4no, respectively. The dash-dotted and dashed lines are 
calculations within the IG model with the same hadron set, again for no and 4no. 
Relaxation time is calculated following Eq. (j58p . The dotted line corresponds the 
the lower AdS/CFT bound r]/s = 



5.4 Collisional viscosity in heavy ion collisions 



Above we have studied the specific viscosities of the hadron matter at dif- 
ferent temperatures and baryon densities. In reality a hot and dense system 
formed in a heavy-ion coUision expands toward a freeze-out state, at which 
the components stop to interact with each other. 

We should specially note that the approximations of the slow hydrodynamic 
expansion are violated at the freeze-out stage, which is assumed to be instan- 
taneous within simple hydrodynamical models. More elaborated approach as- 
suming continuous freeze-out decoupling demands some hybrid of kinetic and 
hydrodynamic description e.g., see PU|^ . Rapid processes out of equilibrium 
may lead to an additional dissipation and particle and entropy productions 
[88|89] . These effects are beyond the scope of the relaxation time approxima- 
tion to the quasiparticle Boltzmann equation used in this work to evaluate 
the collisional viscosity. In addition, there are other sources contributing to 
viscosity, see e.g. Ref. [IH] and also Appendix C. Here we use the phenomeno- 
logical freeze-out curve Tfr(/ibaj,) extracted from the analysis of experimental 
particle ratios in statistical model for many species at the given collision en- 
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1 /2 f 

ergy Sjyjy treating the freeze-out temperature Tfr and chemical potential /ijjg^j, 
as free parameters [86|87] . Thereby extra particle production is phenomeno- 
logically incorporated. We may also hope that including particle production 
effects we also partially include the corresponding effects of the entropy and 
viscosity productions. 

In Fig. [151 C/^ ratios are shown for Au + Au collisions versus the freeze- 
out temperature (which is unambiguously related to the freeze-out chemi- 
cal potential /i^j^j, [87] needed to calculate thermodynamical quantities at the 
freeze-out). As we see, the rj/s ratio decreases monotonously with increase of 
the temperature, being higher than the lower bound 1/4tt but tending to it 
with further increase of Tfj.. The value (/s exhibits a maximum at Tfr ~ 85 
MeV and then tends to zero with subsequent increase of Tfr. As has been 
emphasized above, at T >100 MeV values of the shear and bulk viscosities 
become comparable in value, {jj/ s)^^ ~ 2(C/s)fi.. 
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Fig. 15. Specific shear and bulk viscosities calculated in the SHMC model for central 
Au+Au collisions along the freeze-out curve (at T = Tf^) |87] for the baryon 
enriched system. The dotted line is the lower AdS/CFT bound r//s = I/Att |31j . 

In Fig. [ini the rj/s ratio calculated in our SHMC model (solid line) is plotted 
as a function of the collision energy ^/snn of two Au-|-Au nuclei. The result 
for the IG model with the same hadron set as in SHMC model is plotted by 
the dash-dotted line. We note that for yj/vw ~ 3 the SHMC results prove to 
be very close to the IG based model ones (with the same hadron set as in 
SHMC model), since the freeze-out density is rather small and the decrease 
of the hadron masses occurring in the SHMC model is not important. The 
results for the hadron hard core gas model (the van der Waals excluded volume 
model) [79] at two values of the particle hard core radius r are shown by dashed 
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Fig. 16. The specific shear viscosity calculated for central Au+Au collisions along 
the chemical freeze-out curve j87j within the SHMC model as a function of the 

1/2 

collision energy (solid line). Dashed and short-dashed curves are the results of 
the excluded volume hadron gas model [79] with hard-core radii r =0.3 and r =0.5 
fm, respectively. The dot-dashed line corresponds to the IG model with the same 
set of hadrons as for the SHMC model. 

and short-dashed lines. In all cases for ^/SNN ~ 2 GeV the ratio r]/s decreases 
along the chemical freeze-out line with increasing collision energy and then 
flattens at a/J/vw <;10 GeV, since freeze-out at such high collision energies 
already occurs at almost constant value of Tfr ^ 165 MeV. The shear viscosity 
of the non-relativistic Boltzmann gas of hard-core particles [79] is oc VmT/r^. 
Since Fermi statistical effects are not included within this model, the shear 
viscosity, rj, decreases with decrease of T. Nevertheless the 77/5 ratio increases 
and diverges at low energy/temperature, as the consequence of a more sharp 
decrease of the entropy density compared to rj, see Fig. [161 As follows from 
the figure, the smaller r is, the higher rj/s is in the given excluded volume 
model. For ^/snn ~ 4 and r ^ 0.7 fm the rj/s ratio is expected to be close to 
the values computed in the IG and SHMC models. We also note that in the 
whole range of the considered -^/snn the SHMC results proved to be very close 
to the IG-based model expectation (with the same hadron set as in SHMC 
model) since the freeze-out density is small. 

Recently an interesting attempt has been undertaken in [92] to extract the 
shear viscosity from the 3-fiuid hydrodynamical analysis of the elliptic flow in 
the AGS-SPS energy range. An overestimation of experimental elliptic flow f 2 
values obtained in this model was associated with dissipative effects occurring 
during the expansion and freeze-out stages of participant matter evolution. 
The resulting values oi rj/s vary in interval ?7/s~l — 2in the considered 
domain of ^Js^n ~ 4 — 17 GeV (corresponding to temperatures T 100 — 115 
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MeV). Authors consider their result as an upper bound on the rj/s ratio in 
the given energy range. Note that the mentioned values are much higher than 
those which follow from our estimations given above and presented in Figs. 
[T^ and [ini One should additionally mention that the local equilibrium in 
the three-fluid models is established in every fluid but not among fluids in a 
common cell, whereas other models considered above imply a local thermal 
equilibrium. 



Other microscopic estimate of the share viscosity to the entropy density ratio 
for the relativistic hadron gas based on the Kubo formulae was performed in 
Ref . |93] . Calculations were carried out by simulation of a system evolution in 
a box with appropriare boundary conditions using the UrQMD code, where 
55 baryon species and their antiparticles and 32 meson species were included. 
The full kinetic and chemical equilibrium is achieved at T =130 and 160 MeV, 
respectively. Treating these temperatures as kinetic and chemical freeze-out 
ones and comparing them with those in Fig. [151 see that the extracted 
ratio ?7/s >1 exceeds the SHMC result by a factor of 5. Introducing a non-unit 
fugacity or a finite baryon density allows one to decrease the ratio twice but 
nevertheless it is still too high as compared to both the SHMC result and 
the lower bound 'q/s = l/Air. Analyzing their result authors [93] conclude 
that the dynamics of the evolution of a collision at RHIC is dominated by 
the deconfined phase (exhibiting very low values of r]/s) rather than by the 
hadron phase. Note however that in-medium effects in the hadron phase are 
not included into consideration in the UrQMD model though, namely, these 
effects result in the required decrease of the rj/s ratio in our SHMC model. 



Similar dynamical estimate has been done recently in Ref. [HI] where addition- 
ally to real hadrons the Hagedorn states were considered. Inclusion of heavy 
Hagedorn clusters decreases the rj/s ratio at T <140 MeV as compared to 
the UrQMD result [93] cited above but the obtained value 77/5 ~ 0.3 — 0.4 is 
still higher than the lower AdS / CFT bound and the Hagedorn gas statistical 
estimate for /ibar =0 [5U] (see Fig. |H]). It is worthy to note that for a finite 
chemical potential (or baryon density) the specific shear viscosity in [53] is 
getting lower but in contrast with our SHMC model (see Fig. [T^ the lower 
AdS/CFT bound is not reached. 
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6 Shear and bulk viscosities in the quark phase 



6. 1 Estimates of the relaxation time 



In the framework of the quasiparticle approximation the shear and bulk vis- 
cosities in the QGP phase can be found with the help of Eqs. ( 13^ . These 
equations should be supplemented by expressions for the coUisional relaxation 
time. Like in the hadron system, the coUisional relaxation time r^, a = {q,g}, 
depends on the momentum. To simplify our calculations of the transport coef- 
ficients we will use for r values fa, where fa is estimated with the thermal aver- 
aged cross-sections describing total elastic scattering of medium constituents. 
Unlike the hadronic cross section the quark/gluon (parton) elastic scattering 
cross section is not measurable and should be evaluated in a model. Really, 
it is a specific problem. The in-medium cross sections for quark- ant iquark, 
quark-quark and antiquark-antiquark scattering processes were studied in de- 
tail in Ref. [53] within the NJL model for two different flavors, including l/N^ 
next-to-leading order corrections. These results incorporate dominance of the 
scattering on large angles and take into account a possible occupation of par- 
ticles in the final state, see Ref. [96] . 

The QCD calculations of the relaxation time f of partons already in the lowest 
order in the running coupling constant g require summation of infinitely many 
diagrams. Resummation of the hard thermal loops results in the width f~^ of 
partons ~ g'^Tln{l/g) P^. Thus, the following parametrization was used for 
gluons [271198] 



r,-^ = 2iV.^ln^, (60) 
" Sir g'^ 

and similarly for quarks 



. 1 Nl-lg^T^ 2c 

Ta=2— — In—, (61) 

where the running coupling constant is given by 



(lliV,-2iV;)ln(A(r-r,)/T,)2 ' ^ > 

which permits an enhancement near Tc (by construction, the hadron-quark 
phase transition in this model is of the first order). Parameters of the effective 
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coupling Tg and X/Tc were adjusted to the lattice QCD EoS similar to what 
we did above for the two-phase EoS. One uses T^/Tc =0.46 and the tuning 
parameter c is determined from the condition g^(T = T^) = [27|^ that 
results in the divergence of the relaxation time at (the cross section is zero 
at this point). 
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Fig. 17. Temperature dependence of the relaxation time Tq for /ibar — 

0. The solid 

line shows result of a calculation within the NJL model [96]. The quark relaxation 
time Tq estimated according to Eq. (I6ip at T/Tc > 1 is given by the dash line. For 
comparison the relaxation time Ttttt obtained following Eq. (j4ip . but only for the 
gas of massive free pions at T/Tc < 1, is plotted by the dash-dotted line. 

The temperature dependence of fg for /ibar = within the NJL model [96] is 
shown in Fig. [T7]by the solid line. Both T < and T > T^. regions are covered 
in this model. However the NJL model, in which the medium is assumed to 
be composed only of heavy constituent quarks, can hardly be applied for 
the description of the T < Tc phase. Therefore at T < Tc for this phase we 
suggest to use the developed above hadronic description, where collisional time 
is calculated following Eq. (jUj) with experimental cross-sections for hadronic 
species. In the case /ibar = the main contribution is given by the pion term 
(see dash-dotted curve). In Fig. [17] we also show the relaxation time fq for 
T > Tc following ( l6Ti) . see the dashed line in figure. In order to obtain the 
relaxation time for gluons, Tg, we should multiply fg by the group factor 4/9. 

6.2 Collisional viscosity in the quark- gluon phase 

In the large Nc limit the specific shear viscosity in the QGP phase can be 
inferred from calculations of [35j 
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V S / QGP 



'1 + 3.974A 



69.2 

'ln(26/(^2;Y^(l + 0.5O) 



(63) 



with ^ = Nf/Nc- Obviously this ratio in the quark/gluon phase at finite large 
Nc differs from that calculated for the hadronic phase. Therefore, one may 
expect a jump in the temperature dependence of the ?7/s-ratio at the crossing 
of the phase boundary of the first-order phase transition. In our HQB model 
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Fig. 18. The T-dependence of the shear (left panel) and bulk (right panel) specific 
viscosities within our two-phase SHMC-HQB model (solid lines) for ^bar =0. The 
results for the original NJL model [96] are plotted by dashed lines. The rj/s and C/s 
calculated following Eqs. (j52|) . (|53|) using quark masses from the HQB model and 
the relaxation time from (|6ip are shown by the dash-dotted lines. The lower bound 
for the reduced shear viscosity, l/47r, is given by the horizontal dotted line. 

strange quarks and gluons are very massive. Therefore their contribution to 
transport coefficients and entropy density is minor. Main contributions arise 
from u- and d-quarks. 

In Fig. [iHl rj/s and (/s for hadron and quark-gluon phases are shown in a 
broad temperature range for /Xbar = 0. Solid lines present results of the two- 
phase SHMC-HQB model. Below Tc the ratios are the same as shown above in 
Fig. IHl Calculations done in the NJL model [96] are shown by dashed lines. In 
both cases the quark collision relaxation time is that shown in Fig. [T7] by the 
solid line. We see that the original NJL model [HS] gives a continuous smooth 
line (a crossover) for /ibar = 0. In our SHMC-HQB model there is a jump at Tc 
in both the rj/s and (/s ratios. This jump is a particular property of the first 
order phase transition. For /ibar =0 in our two-phase SHMC-HQB model the 
rj/s ratio increases with T more sharply after the jump has occurred compared 
to that in the original NJL model (see Fig. ITSl) . 



Dash-dotted curves in Fig. [TS] demonstrate that the specific viscosities diverge 
for T — 7- Tc + since f — ?■ oo provided the nonperturbative estimate of the 
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relaxation time ( 16T|) is used. With increase of the temperature above Tc the 
temperature dependence of these ratios flattens. The gluon rj/s and (/s ratios 
exhibit in this similar divergent behavior near T^, since the relaxation 

time (pi also diverges at + 0. As follows from Eq. (jSSD, in the QGP phase 
rj/s oc where in accordance with fl62l) the coupling constant is expressed 
via g"^ (X 1/ ln(T/AT), At is proportional to the scale parameter of QCD. Sim- 
ilar minimum close to Tc from the QGP side arises for the rj/s for the gluon 
quasi-particle excitations in a phenomenological quasiparticle model |100] in 
agreement with lattice QCD data in gluodynamics [5^4U|41j . Thus we con- 
clude that both the magnitude of the jump and the rj/s behavior for T > T^. 
are essentially model-dependent. Generally, similar conclusion can be done 
from comparison of different curves for (/s in the right panel of Fig. [181 

Note that for the hadron phase the NJL model noticeably overshoots the 
SHMC model results for both rj/s and (/s specific viscosities. However one 
should be careful comparing the NJL and SHMC-HQB results. The NJL model 
treats the chiral phase transition while in the SHMC-HQB model we deal with 
deconfinement. Thus, critical temperatures have different meaning in these 
models. Their values are also different (but the latter problem can be avoided 
provided one uses T/Tc as an argument). Moreover in the NJL model one 
continues to deal only with u and d quark degrees of freedom at T < 
whereas the SHMC model involves many hadron degrees of freedom whose 
masses essentially decrease at T ~ Tc. Thereby, the entropy density in the 
SHMC model is essentially higher and the specific viscosity is accordingly 
lower than those in the NJL model [96]. The shear viscosity-to-entropy ratio 
in the latter model is much higher than its lower bound l/47r at T T^. 

One should note that in our model the deconfined partonic phase is described 
in a simplified way as an ideal massive gas but relaxation times were taken from 
the models including interaction. A consistent approach requires an account 
for interactions between the constituents and the extracted effective couplings 
should enter the estimate for the transport cross section. A step toward this 
direction was undertaken in Ref. |1U1] within the generalized classical virial 
expansion formalism. The corrections to a single particle partition function 
were calculated starting from a parton interaction potential whose parameters 
are fixed by lattice thermodynamical quantities. The dependence of the r]/s 
ratio on the temperature dependence of the strong interaction coupling has 
been studied. Results prove to be sensitive to the choice of the interaction. 

There exist many other model estimates of the rj/s ratio for T > Tc for yU-bar =0. 
As noted in the Introduction, AdS/CFT [2T] predicts rj/s = l/(47r). Lattice 
Monte Carlo calculations are available only for gluodynamics. For the SU(3) 
gauge theory a robust upper bound rj/s < 1 was obtained with the best es- 
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timate rj/s = 0.134 at T = 1.65Tc [H]. A slightly larger value of the shear 
viscosity was derived in the lattice QCD simulation in the quench approxima- 
tion for a hot gluon gas [SH]: ?7/s — O.l-i-0.4 for temperature in the region 
1.4 < Tc < 1.8. At higher temperatures 7] increases by two-three orders of 
magnitude. 
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Fig. 19. The T-dependence of the shear (left panel) and bulk (right panel) specific 
viscosities within our two-phase SHMC-HQB model (solid lines) for ^i^^^ =990 
MeV, corresponding to the critical end point. The NJL model results j96j for ^^a^^ 
and for //bar slightly above ^^^^ are plotted by long-dash and short-dash lines, 
respectively. 

Model predictions for the rj/s and C,/s ratios for systems with /ibar 7^ are 
presented in FigHH The NJL model [96] at the critical end point (/ibar =990 
MeV) shows a smooth curve (long-dash line) with a sharp bend at for rj/ s 
and C/s. Certainly, our two-phase model again exhibits a jump at (see solid 
lines). Similar jump will appear in the NJL model if //bar is larger than 990 
MeV, i.e. when the system reaches the region of the first order phase transition 
(see short-dashed lines). 

Here the following remark is in order: As follows from all above estimates, the 
bulk viscosity is smaller (but not much smaller in the hadron phase at finite 
baryon density) than the shear one in the whole temperature interval, may be 
except the vicinity of the critical point. In the latter region soft modes may 
additionally contribute significantly increasing the bulk viscosity. 

In spite of a variety of existing models, none of them is reliable for the de- 
scription of the vicinity of the critical point. In this respect, it is highly de- 
sirable to get experimental information about viscosity realized in heavy ion 
collisions to compare it with the model predictions. The shear viscosity is in- 
timately related to the observed elliptic flow V2 and viscous hydrodynamics 
provides a tool to extract rj/s from experimental data. However, calculation 
uncertainties are significantly affected by details of the initialization of the 
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initial hydrodynamic state, behavior of the bulk viscosity and the speed of 
sound near the quark-hadron phase transition point, chemical composition 
and strong non-equilibrium effects during the late kinetic stage of the evolu- 
tion of the hadronic state. From the study of transverse momentum fluctu- 
ations at the RHIC energies one gets 77/3 ~ 0.08 ^ 0.3 in Ref. |102] . Close 
estimates rj/s ~ 0.08 -r- 0.3 and rj/s ~ 0.09 ^ 0.15 were obtained from V2 
analysis in Refs. |103j and |104] . respectively. In its turn, studies of the heavy 
quark energy loss and meson production yield ///s ~ 0.1 -r- 0.16 |105j and 
?7/s ~ 0.07 ± 0.003 ± 0.014 |106] . respectively. Recent comprehensive analysis 
of the hydrodjTiamic simulations |107] has established an upper limit 



QGP < 5 X ^ ' ^^^^ 



to be rather close to the low AdS / CFT bound. Note that the presence of C 7^ 
can also result in a suppression of the elliptic flow V2 competing at T near T^. 
with shear viscous effects |108] . Thus, the above constraint on 77/5 can be even 
stronger. 

Finally, one should emphasize that practically in all viscous hydrodynamic 
calculations the specific ratio 77/5 is considered as a time-independent quan- 
tity, that is not the C3iS6, clS follows from the above-presented dependence of 
kinetic coefficients on temperature and baryon density. In addition, it is tac- 
itly neglected by retardation effects in the formation of irreversible current. 
Some exceptions are estimates |109fll0j . In Ref. |109j viscosity effects were 
studied with two different T-dependent shear and bulk viscosities (treated as 
free parameters) in the QGP and hadronic phases. The values of the used 
transport coefficients noticeably differ from those presented in FigsJT8]and[T9l 
Following |109j the elliptic flow coefficient is significantly reduced due to vis- 
cosity effects both in the plasma and in the hadron matter. A more realistic 
parameterization of the temperature dependence of shear and bulk viscosity 
coefficients was considered in Ref. [110] and additionally relaxation times were 
included in order to explain elliptic flow at RHIC It was argued that close 
to ideal behavior observed at RHIC energies may be related to a non-trivial 
temperature dependence of these transport coefficients. So, extraction of the 
shear viscosity from analysis of the initial QGP phase is a difficult problem 
which can be reliably addressed only after precise constraining the freeze-out 
conditions. Thus, at present there exist many different estimates of viscosity 
coefficients, being not consistent with each other and only indirectly related 
to the observed flow effects. Thus no deflnite conclusion can yet be done about 
temperature and density dependence of the transport coefficients in the QGP 
and hadron phases. 
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7 Beyond quasiparticle approximation 



7. 1 Generalization to finite particle mass-widths 



Expressions (HHl) . ( I49l) for tlie shear and bulk viscosities were derived within 
the quasiparticle approximation, i.e., the particle mass-width F = — 2ImS^ 
was put zero in the retarded Green function G^. Here is the retarded self- 
energy. The spectral density [A = — 2ImG'^) is then expressed through the 
5-f unction. 

Let us start with the consideration of a fermion spin 1/2 resonance, "f. The 
spectral function satisfies the sum rule: 



in. 

4 JO 



7o 
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{+){Po,p} + Al^)iPo, 



-V) 



dpo 
2tt 



(65) 



where 70 is the corresponding Dirac matrix, and subscripts (±) specify par- 
ticle and antiparticle terms. The trace is taken over spin degrees of freedom. 
Simplifying the spin structure we introduce the spectral density 



A' = ^Tr{ 70^;+)}. (66) 

Consider the dilute matter assuming that the spectral function depends only 
on the variable s = Pq —p^ > 0, see [111] . To do the problem tractable, instead 
of solving a complete set of the Dyson equations, we will use a simplified 
phenomenological expression [55] : 



,f e 2po [T^js) + 6] 

A = A 2po = o ^ , c = const, 67 

^ {s-mr y + [T'{s) + 6]y4 ^ ' 

with 6 —7- +0 introduced to easier perform the quasiparticle limit. The value 
r^(s) is the width. Phenomenological expression for its s-dependence is pre- 
sented in Appendix B. 

For charged bosons the spectral function has the form0 



^ As before, by the charge we mean any conserved quantity like electric charge, 
strangeness, etc. 
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e [rbos(s) + ^] 



^bos — ^ part*N2N2 ■ , no/^' v""^' 



Replacing Fbosl*) with T{s) and Abos with A, we may use Eq. (!82|) of Appendix 
B as a phenomenological parameterization of the width. The quantity ^ is 
introduced to fulfill the sum-rule in the form 



We continue to study dilute matter assuming that the spectral function de- 
pends only on the s- variable. Only in this case one may consider a single spec- 
tral function for vector mesons, like u and p, whereas in the general case one 
should introduce transversal and longitudinal components. Now we conjecture 
expression for the shear viscosity which generalizes quasiparticle expression 
obtained above to the case of the finite mass-width 



Here we have put = 2pQ/Ta{s), Ea = \/(mr'^**)2 +fp^ 



Derivation of this expression will be given elsewhere. In the nonrelativistic 
approximation Eq. (ITOl) coincides with the corresponding expression derived 
in [TT2] . 

In the quasiparticle approximation 



^^27r5(.-(mrT), ^ 2^6(3 - {mlT?), (72) 

and we return to Eqs. ( HHl) - ( l49l) with = 2Ea/Ta{Ea). Here the relaxation 
time for the a-quasiparticle is expressed through the width, rather than the 
cross-section. Bearing in mind the optical theorem, these two quantities should 
coincide within the quasiparticle approximation. This statement however re- 
quires an additional check that we do not perform in the given work, since our 
aim here is just to estimate possible effects of finite widths on the viscosities. 
Moreover in the case of a finite width the off-mass shell Ta{po,p) enter our 
expressions. 
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Fig. 20. Temperature dependence of the ratio of the shear viscosity for the resonance 
calculated in the constant width approximation at different values of f = Fomres to 
that computed in the quasiparticle approximation. In both (resonance and quasi- 
particle) cases we use the same (for the given Fq) energy-momentum independent 
values of the relaxation times and m-res = rn-A- Lines show different values of the 
width. 

In Fig. |20l we show the temperature dependence of the ratio of the shear 
viscosity calculated for the resonance (mres = ""^a =1.232 GeV) in the con- 
stant width approximation (F = Fomj-es) with the constant relaxation time to 
that computed in the quasiparticle approximation with the same relaxation 
time. Realistic value of the A-resonance width in the resonance maximum is 
-pn.rei. ^ ]^20 MeV that corresponds to the choice Fq = 0.24 GeV. In order to 
understand the dependence of the shear viscosity on the value of the width 
we vary the value Fq. We see that width effects are quite important for the 
description of the resonance characteristics at low temperatures T < 50 ^ 100 
MeV, whereas at higher temperatures one may use quasiparticle approxima- 
tion. At low temperatures contribution of the resonances to the transport 
coefficients and to the thermodynamical quantities is suppressed compared to 
the nucleon and pion contributions. However, since even for very small width 
(for Fq = 0.01 GeV corresponding to F'^''"''^' ~ 5 MeV) the quasiparticle ap- 
proximation for calculation of the viscosity fails for T < 50 MeV, and since 
nucleons have some width, we may conclude that calculations of the transport 
coefficients at low temperatures should be performed with taking into account 
width effects. 

It was argued in |113] that in a theory, where the lowest energy excitations (not 
quasiparticles) form a continuum, the rj/s ratio has no lower bound in contrast 
with predictions of the AdS/CFT correspondence. To check this statement we 
calculated rj/s ratios. For the free resonance gas the entropy density is given 
by [TTTITTl], 
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Fig. 21. Temperature dependence of the ratio of the entropy densities for the 
resonance calculated with constant width and in the quasiparticle approximation, 
TTT-res = ru/^. The curves are drown for different values of the width. 
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Fig. 22. Temperature dependence of the specific shear viscosity for the A resonance 
calculated in the constant width approximation, m^^s = TnA- The quasiparticle 
result is plotted by dots. 

where = t{1 T F;^"^) ln(l T F^"^) - F^'^lnF^'^ and F^'i is determined by Eq. 

In Fig. [21] we present the ratio of the entropy density computed with taking 
into account width effects to that calculated in the quasiparticle approxima- 
tion. As for the viscosity, in case of small width (see solid curve for Tq = 0.01 
GeV) the ratio s^'^^/s'^^ deviates from unity only for low temperatures (at 
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T < 50 MeV). For a larger width the ratio begins to significantly differ from 
the unity in the whole temperature interval. 

In Fig. [22]the temperature dependence of the -q^'^^ / s'^'^^ ratio is shown for mres = 
mA and for different values of Fq. We see that the dependence on the value 
of the width is rather moderate. The ratio calculated in the quasiparticle 
approximation behaves in a similar way but the maximum position is slightly 
shifted toward a higher temperature. For large temperatures the quasiparticle 
limit is achieved and, thereby, within our model the rj/s ratio reaches the limit 
l/47r at T ~ ~ 180 MeV. 

7.2 Mean- field bulk viscosity term 

The sources of the shear and bulk viscosities considered above are associated 
with collisions between in-medium excitations (quasiparticles and particles 
with widths). Another source of the bulk viscosity may come from a dissipation 
of the soft modes related to a slow dynamics of the mean fields provided the 
system is located in the vicinity of the critical point of a second order phase 
transition or at a weak first order phase transition. Note that in our SHMC 
model the soft modes relate to the hadron masses decreasing towards the 
critical point T^. Such a contribution was first introduced in [ITPS] and then 
discussed in [1^50f51] . Within the quasiparticle approximation used above, 
the mean fields do not dissipate. Thereby, there are no linear in terms in 
P3l) responsible for an additional bulk viscosity term. In reality there may 
exist such a dissipation. We demonstrate possible effect of a dissipation on the 
bulk viscosity in Appendix C. However we should note that considering the 
soft mode viscosity regime one assumes that the time characterizing the soft 
mode evolution is sufficiently long. Estimations [51] show that soft modes 
may not develop during the time, at which the system trajectory in the course 
of the heavy ion collision passes the vicinity of the critical point. Thus the 
particle collisions look a more adequate source of the viscosity in heavy ion 
collisions than a source due to the soft mode evolution. Therefore in the given 
paper we focused on discussion of the collisional source for the viscosities. 



8 Conclusions 

In this paper, we analyzed different approximations for the calculation of the 
shear and the bulk viscosities for the hadron and the quark-gluon phases. 
General expressions for shear and bulk viscosities are derived in the relaxation- 
time approximation for a system described by the quasipartcle relativistic 
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mean-field theory with the scaling of hadron masses and couplings (SHMC). 
The EoS of the SHMC model fairly well reproduces global properties of hot 
and dense hadron matter including the temperature region near and the 
lattice data for < T < 1.3Tc, provided all coupling constants g^rb are strongly 
suppressed except for the nucleons. 

We compared kinetic coefficients calculated in the SHMC model with those 
calculated in other models of the hadron phase. At /Xbar = a general fall 
off of the specific shear viscosity with temperature in hadronic phase is a 
common feature for all models but quantitative values are somewhat different 
particularly showing different behavior in the temperature region near the 
transition temperature Tg. The rj/s ratio in the SHMC model is closer to that 
of the excluded volume hadron-resonance gas model with the hard-core radius 
r ^ 0.7 fm [79] and practically coincides in the near-critical region with a 
resonance gas model result that includes Hagedorn states |115j . 

For /ibar 7^ transport coefficients qualitatively behave similarly to the case of 
baryonless matter reaching the lower AdS/CFT bound for the specific shear 
viscosity rj/s = I/Att with increase of T. For the higher baryon density this 
limit is reached at smaller T. In particular, with increasing freeze-out tempera- 
ture Tfr (for central Au-|-Au collisions), the rj/s ratio undergoes a monotonous 
decrease approaching values close to the AdS/CFT bound at T ~ TcQ while 
the (/s ratio exhibits a maximum at Tfr ~85 MeV. In a broad temperature in- 
terval the rj/s and (/s ratios are not small and viscous effects can be noticeable. 
The viscosity values at the freeze-out can be transformed into dependence on 
the colliding energy y^J/vw (for central Au-|-Au collisions). When the collision 
energy decreases, the rj/s goes up. The high-energy flattening of the a/J/vtv 
dependence occurs at quite low ri/s < 0.2. It implies that a small value of rj/s 
required for explaining a large elliptic flow observed at RHIC could be reached 
in the hadronic phase. This might be an important observation which we have 
demonstrated within the SHMC and the ideal gas (IG) models with the same 
hadron set. The rj/s and (/s ratios at the freeze-out curve calculated within 
the SHMC and IG models agree with each other; however, rough estimates 
P^93|94j oi rj/s extracted from comparison of dynamical model calculations 
with experimental data overshoot these theoretical expectations by a factor 
of 2-5. Note that rapid processes at the freeze out may result in an additional 
increase of the viscosities and the entropy production. These effects are not 
incorporated directly in the present consideration but the use of experimen- 
tal freeze-out temperature and chemical potential allows one to hope that 
indirectly they are effectively included. 

We remind that there exist rather general arguments that rj/s should have a 
minimum at the QCD phase transition critical point similar to that exists for hehum, 
nitrogen and water, see |37j . 
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Using the Gibbs conditions we extended our approach to higher temperatures 
combining the SHMC model description of the hadron phase with that of the 
heavy quark bag (HQB) model for the quark-gluon phase. It was demonstrated 
that this two-phase SHMC-HQB model is in reasonable agreement with QCD 
lattice data at /ibar =0 and /ibar 7^0 at all temperatures except a vicinity of the 
critical point. The finite-size and nonequilibrium effects of the phase mixture, 
possibility of a rapid phase transition without a phase mixture may lead to 
significant additional dissipation provided the system trajectory passes the 
vicinity of the critical point. Moreover one might be need to improve the model 
in order to include possibility of the crossover for yUbar = 0. Not considered in 
our work these effects require a special careful analysis. 

In the two-phase SHMC-HQB model with the first order phase transition there 
appears a jump in both 77/5 and (/s at the critical temperature at any /Xbar- 
For T > Tc at yUbar = 0, 77/5 grows with increase of the temperature, whereas 
C/s approximately stays constant. The magnitude of both ratios differs sig- 
nificantly from those estimated in the NJL model [59|96] . Comparison of our 
approach with the NJL model [5^1^ demonstrates that behavior of transport 
coefficients near the critical point depends on the type of the phase transition. 

It is worthy to note that the V2 analysis [116j indicates different values of 
•q/s for peripheral and central collisions. Therefore, it would be interesting to 
perform hydrodynamical calculations using the T — /ibar dependent transport 
coefficients rather than constant ones. Necessity of such an approach was re- 
cently emphasized in [116] . In subsequent works we plan to use the SHMC 
model EoS with the derived transport coefficients for this purpose. 

We also estimated influence of the finite width effects on the transport coeffi- 
cients. For this aim we evaluated the shear viscosity and the entropy density 
at non-zero width. Our findings show that width effects are important for low 
temperatures (for T < (50 -7- 100) MeV), whereas for higher temperatures one 
may use the quasiparticle approximation. 

As was demonstrated in [501I5T] the (/s ratio may have a sharp maximum at 
the critical point. This statement is based on the assumption of the existence 
of soft collective modes at temperatures near T^. The enhancement of the ratio 
(/s near the critical point could have interesting consequences [T^17II43|45|117] 
provided the system evolves very slowly. However it is unlikely that in the 
course of relativistic heavy ion collisions the system spends enough time in 
the vicinity of the critical point to develop slow soft collective modes, see 
[M|118j . Moreover effective hadron masses responsible for the softness of the 
modes, being calculated within our SHMC model, although decrease toward 
the critical point do not drop to zero. Thus the quasiparticle-like coUisional 
viscosity estimates look more relevant than those follow from the considera- 
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tion of the slow soft mode dynamics. One should remind once more that other 
sources for an increase in the bulk viscosity may exist, as discussed in |19]. 
Therefore further studies of the given problem are still required. 
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Appendix A. Necessary formulas for derivation of transport coeffi- 
cients. 

Here to take derivatives in ( 142|) we follow the line of Refs. |15|59] . The only 
difference is that the strangeness conservation is additionally incorporated. 

The energy density e and the baryon and strangeness charge density (nbar and 
ristr) conservations can be expressed as 



de 
dt 

'9nbar 

dt 
dn^tr 
dt 



V ■ 



dP dP ^ 

+ fJ'strT; V ■ n, 

dfistrj 



-^barV ■ U 
-ristrV ■ u -- 



dP 



dfihar 

dP 



V ■ u 



(74) 
(75) 
(76) 



Here we set m = keeping only derivative terms of u. The pressure is expressed 
in terms of T, /i variables. Further expressing P in terms of e, n variables and 
using Eqs. (!74|) - (176|l we find 



dP[e,n] _ dPde 
di ~ 'd^di 



dP Sribar dP dn. 



str 



Sribar dt 

dP fdP 

-TT- T— + /_ibar 



de V dT 
dP dP 



dfihar 

dP dP 



dUstr dt 

dP dP 

"I" Z^str" 



str , 



^ribar ^yUbar dn^tr djJ, St 



V ■ u. 



(77) 



On the other hand 
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dP[T,fi] _dPdT , dP a/ibar ^ dP a/istr 



dt 

and thus 



dT dt 9/ibar dUstr dt 



9t [ de . 



V -M, 



dt 



'dP] 
de 



9P ~ 



dn. 



str , 



V •£ 



(78) 

(79) 
(80) 
(81) 



We use these equations to construct variation of the energy- momentum tensor 
5T'^ in Eq. (gS]). 



Appendix B. Phenomenological expression for the width of the res- 
onance. 

The phenomenological expression for the s-dependence of the width T^{s) can 
be easily recovered in the near threshold region, see 




0{s - sth) 



12) 



Here a = 1/2 for the s-wave resonance and 3/2 for the p-wave resonance, is 
the fermion resonance mass in vacuum, and Sth is the threshold value of s; Fq is 
the constant of the dimensionality of the energy. In order to use this expression 
outside the near threshold region one introduces the form-factor -F(s), the sq 
is the cut-off constant and the power /3 > 1 -|- a/2. The parameters can be 
adjusted to satisfy experimental data. The factor ^ is introduced to fulfill the 
sum-rule: 



°° ds -f 

— Ai = 1 . 
27r 



(83) 



For charged bosons the spectral function follows the sum-rule, 



r°° ds 
lo 2^' 



^bos — 1 • 
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Appendix C. Contribution of soft collective modes to the bulk vis- 
cosity 

As follows from the particle data, in vacuum the width of the a meson is 
very large, Fg- ~ 400 MeV. In hot baryonic matter the width can additionally 
increase due to a collisional broadening. Supposing a = acq+Sa in the equation 
of motion for the mean field, where acq is the equilibrium value of the mean 
field, we arrive at the equation for the fiuctuating mean field part 6a: 



im*^'[T,^^]r6a = -T^^, (85) 

l/Fo- is the relaxation time for the given process. We suppressed the second 
space-time derivative terms but included the first time-derivative dissipation 
term in this equation with a pre-factor Fg-. From Fig. [1] we see that for the zero 
baryon chemical potential in the SHMC model with suppressed couplings gab, 
except for the nucleons, which we use in the given paper, the effective mass 
gnificantly decreases at r ~ T^,. 

Using smallness of the time-derivative term in fl85p we find 



where cxioceq is the mean field value at the local equilibrium. On the other 
hand, 

= , -77 = Sloc.cq V ■ U, (87) 

dt ds dt ds 

where we used approximate conservation of the entropy (up to terms propor- 
tional to transport coefficients assumed to be small). Since the bulk viscosity 
contribution of the mean field is related to the pressure as Cmf V ■ m = —5Pmfi 
we find 



/9Pmf\ 6a (dPMY\ \ r da A 

and [^jioc.cq = [%l ^]ioc.oq IS takcu at the condition that equations of motion 
(fTT!) are fulfilled. 

There exist arguments that the p-meson also becomes a broad resonance at T 
near Tc, see [119] , and its effective mass decreases with increase of T towards 
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Tc- The p-mean field appears in isotopically asymmetric nuclear matter. In 
this case, there may appear an extra contribution to the bulk viscosity. A 
similar contribution from the fluctuation of the loq field is probably suppressed 
compared to that from the cr, since F^j <^T„. 
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